Critical exponents of four-dimensional random-field Ising systems 
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The ferromagnet-to-paramagnet transition of the four-dimensional random-field Ising model with 
Gaussian distribution of the random fields is studied. Exact ground states of systems with sizes 
up to 32* are obtained using graph theoretical algorithms. The magnetization, the disconnected 
susceptibility, the susceptibility and a specific heat-like quantity are calculated. Using finite-size 
scaling techniques, the corresponding critical exponents are obtained: /3 = 0.15(6), 7 = 3.12(10), 
7 = 1.57(10) and a = Q (logarithmic divergence). Furthermore, values for the critical randomness 
tic = 4.18(1) and the correlation- length exponent v — 0.78(10) were found. These independently 
obtained exponents are compatible with all (hyper-) scaling relations and support the two-exponent 
scenario (7 — 27). 
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I. INTRODUCTION 

Phase transitions of pure systemj^'l are already rela- 
tively well understood. The critical behavior of all phys- 
ical quantities can be described via critical exponents. 
These exponents are related through (hyper-) scaling re- 
lations to each other, so that only two independent expo- 
nents remain. In contrast, phase transitions in systems 
with (quenched) disordem exhibit many puzzles and are 
still far from being understood. 

In theoretical physics, the random- field Ising magnet 
(RFIM) is a wideWj-studied prototypical disordered sys- 
tem. It is believedO'El to be in the same universality class 
as the diluted antiferromagnet in a field, which can be 
studied experimentallyQ. 

For a while it was thoughtLJElEI that the critical behav- 
ior of the d-dimensional RFIM is equal to that of the 
d—2 pure ferromagnet. This would imply that the d = 3 
RFIM exhibits no ordered phase. This is not true, as 
has been showH|4ftter rigorousljtj. In the meantime, a 
scaling theorylliUla'tj for the RFIM was developed, where 
the dimension d has been replaced by d — in the hyper- 
scaling relations, 9 (sometimes called also y) being a ttiird 
independent exponent.. iii contrast to the pure case. An 
alternative approacnlJlia leads to the consequence that 
9 is not independent but related to the exponent 77 de- 
scribing the divergence of the (disconnected) susceptibil- 
ity via 9 = 2 — rj. Further evidence for the existence 
of only two independent expoaent was recently found by 
high-temperature expansional3. This was confirmed in 
three dimensions by Monte Ca*kj|-simulationsEZl and by 
exact ground-state calculationstS'ta. But the exponents 
found in these works do not fulfill the scaling relation 
a -t- 2/3 -f 7 = 2, being a,/3,7 = i^(2 — 77) the critical 
exponents for the specific heat, the magnetization and 
the susceptibility, respectively. On the other hand, the 
results obtained in a different way in the most thorough 
ground-state study in d = 3 so farE3 indeed do not vio- 
late this scaling relation. Also,pa-modificd scaling relation 
a -t- 2/3 -I- 7 = 1 was proposed^, but the exponents ob- 
tained in Refs. |l7|Jl^,n9 do not match the new relation 



either. Hence, to obtain more insight into the scaling 
behavior of the RFIM and to improve the knowledge of 
the critical behavior of random systems, here the four- 
dimensional model is studied. 

The RFIM Hamiltonian is given by 



n = -jY, s^SJ - ^(/i,; + H)s, 



(1) 



where the Si — ±1 are Ising spins, J is the interaction 
energy between nearest neighbors, tii = tiei is the ran- 
dom field and H an uniform external field. Here the case 
-ff = is studied, but small values of the external field 
are used to determine the susceptibility. The values 
are independently distributed according a Gaussian dis- 
tribution with mean and standard deviation 1, i.e. the 
probability distribution is 
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Hence, tii — liti is Gaussian distributed with standard 
deviation li. We shall consider finite-dimensional lattices 
with periodic boundary condition and N = L'^ spins. 
The results presented below are for d = 4. 

The canonical phase diagram in zero external field 
[H = 0) of the RFIM in higher than two dimensions 
is shown in Fig. |l|. For low temperatures and small ran- 
domness, the ferromagnetic couplings dominate, hence 
the system exhibits a long-range order. For higher tem- 
peratures, where the entropy dominates, or for higher 
random fields, where the spins are predominately aligned 
parallel to its local random fields, the system is param- 
agnetic. 

For random-field distributions which exhibit a maxi- 
mum at ti — 0, such, as present in the Gaussian case, 
in mean-field thcorycj the phase transition is second or- 
der along the whole phase boundary. Furthermore, the 
rcnormalization group flow along the phas(i-|boundary ap- 
proaches the h = he, T = fixed pointQ, i.e. it con- 
trols the critical behavior. Therefore it can be expected 
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FIG. 1: A sketch of the phase boundary of the random field 
Ising model. The ferromagnetic phase is denoted by "F" and 
the paramagnetic phase by "P" . The critical value of the ran- 
dom field at T = is denoted by he . The lines with arrows at 
both ends indicate the path followed by varying J for some 
fixed value of h and T. 



that the critical behavior along the whole phase bound- 
ary is equal to that at h = hc,T = 0, and all critical 
exponents can be obtained from the zero temperature 
behavior. Working at T = has the advantage that 
exact ground states of the RFIM can be calculated us- 
ing pandern graph theoretical algorithms in polynomial 
timec3'Ot3E3. This avoids equilibration problems often 
encountered in Monte Carlo simulations and, even better, 
allows to study much larger system sizes than before. 



The r = bebavior of the four-dimensional RFIM 
has been studiedEII using exact ground states up to size 
i = 10. In addition to the random critical random field 
he, the exponents v describing the divergence of the cor- 
relation length ^ and (3 for the magnetization were ob- 
tained. In this work, not only much larger systems sizes 
up to = 32'* are considerecL-|but also the critical ex- 
ponents a for the specific heal£j, 7 for the susceptibility 
and 7 for the disconnected susceptibility are obtained in- 
dependently, i.e. without using scaling relations. This in 
turn allows to test (hyper-) scaling relations and to in- 
vestigate the assumptions of two or three independent 
critical exponents. The main results of this paper are 
that the specific heat diverges logarithmically or maybe 
faster and that the two independent exponent scenario is 
supported. 

The rest of the paper is organized as follows. Next, 
the algorithms used to calculated exact ground states are 
briefly explained. In the main section all observables and 
results are presented. In the final section the scaling re- 
lations are checked, the results discussed and a summary 
given. 



II. ALGORITHMS 

To calculate the exact ground ^ 
domnQs^rArapd field H, algorithmpP"^' 




Iljulating the 
and finally 



given ran- 
from graph 

theoryE9Eji23 were applied. To irmlement them, some 
algorithms from the LEDA libraryB were utilized. 

Here the methods are just outlined. More details can 
be found in the literature cited below or in the peda- 
gogical presentation in Ref. |2^. For each realization of 
the disorder, given by the external field H and the val- 
ues {hi} of the random fields, the calculation works by 
transforming the system into a netasint 
maximum flow in polynomial tim(£: 
obtaining the spin configuration {Si\ from the values of 
the maximum flow in the network. The running time 
of the latest maximum-flow methods has a peak near the 
phase transition and divergesHo there like 0{L'^^^) The 
first results of applying these algorithms to random-field 
systems can be found in Ref. ^ In Ref. |l^ these meth- 
ods were applied to obtain the exponents for the mag- 
netization, the disconnected susceptibility and the cor- 
relation length from ground-state calculations up to size 
L = 80. The most thorough study of the ground states of 
the 3d RFIM so far is presented in Ref. |2^. Other exact 
ground-state calculation of the 3d model can be found 
in Refs. |2^,^,^ These techniques hav©jalso already 
applied to small four-dimensional systemg2J. 

Since the algorithms work only with integer values for 
all parameters, a value of J = 10000 was chosen here, 
and all values were rounded to its nearest integer value. 
This discreteness is sufficient, as shown in Ref. |2^. All 
results are quoted relative to J (or assuming J = 1). 

Note that in cases where the ground-state is 
degenerated it is [possible to calculate all the ground- 
states in one sweepllj, see also Refs. For the RFIM 
with a Gaussian distribution of fields, the ground state is 
non-degenerate, except for a two-fold degeneracy at cer- 
tain values of the randomness, where the ground state 
changes, so it is sufficient to calculate just one ground 
state. 



III. RESULTS 

In this work, exact ground states for system sizes L — 
4 to L = 32 for different values of the randomness h 
and with 4 different values H = 0, Hl,2Hl,4:Hl (only 
H = ior L = 32) were calculated. Near the transitions, 
an average over the disorder with the number iVgamp of 
samples ranging from 3200 (L = 32) up to 40000 (L = 4) 
was performed, less samples were used away from the 
critical point, since the fluctuations are small outside the 
critical region. For details, see Tab. |. 

We first concentrate on the case H ^ Q. The simula- 
tions with H > were done to obtain the susceptibility, 
see below. 

As, already mentioned, when studying one single fi- 
nite sample of the disorder {hei} as a function of h, 
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L 




Hl 


4 


40000 


0.025 


6 


20000 


0.02 


8 


7100 


0.012 


12 


8500 


0.006 


16 


4000 


0.003 


24 


10000 


0.0015 


32 


3200 





parameterl 



TABLE I: The maximum number of samples A'samp used, and 
sizes of smallest non-zero uniform field Hl, for each system 
size L. As discussed in the text, the number of samples used 
was larger in the vicinity of the peaks in the susceptibility 
and specific heat than elsewhere. 
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FIG. 2: Average magnetization as a function of random-field 
strength h. For clarity of the plot, only L = 4,8, 16, 32 are 
shown. Error bars (shown for L = 32) are much smaller than 
symbol sizes. Lines are guides to the eyes only. 



the ground state changes only at certain discrete values 
/j(2) _ /jM({ei})^ Hence, quantities like the mag- 
netization are stepwise constant as a function of h. This 
discreteness vanishes, when averaging over the disorder. 
In Fig. ^ the average magnetization per spin 
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is shown as a function of the randomness h for system 
sizes L = 4,8,16,32. The average over the disorder is 
denoted by [. . which is carried out (approximately) 
by repeating the calculation for iVsamp independent re- 
alizations (samples) of the random fields {hei}. As ex- 
pected, for small randomness the ground state is ferro- 
magnetically ordered and disordered for large values of h. 
The curves become steeper with increasing h, indicating 
a phase transition near h = 4.2. 

To study the transition more detailed, the Binder 



(4) 



is calculated, (. . .) being the thermal disorder, (which is 
trivial at T = if the ground state is non degenerate). 
The idea behind the definition of this quantity is that 
in the thermodynamic limit, the distribution of the or- 
der parameter should converge towards a delta function 
(with g{L, h) = 1) in the ordered phase and to a Gaus- 
sian distribution {g{L, h — 0)) in the disordered phase. 
The scaling theory for 2nd order phase transitions as- 
sumes, that the finite-size behavior is governed by the 
ratio L/£^, ^ being the correlation length. At the critical 
point, where the correlation length is infinite, the param- 
eters for different system sizes should intersect, since L/^ 
is sero for all sizes L. In Fig. I the resuh for the 4d RFIM 
is shown, please note the enlarged scale. 




FIG. 3: Binder parameter g{L, h) as a function of the random- 
ness h for different system sizes L. For clarity of the plot, only 
L = 6 is omitted. Error bars are smaller than symbol sizes. 
Lines are guides to the eyes only. 

Indeed all curves intersect near h — 4.2, compatible 
with the result for the magnetization. A systematic shift 
can be observed, which decreases with growing system 
size, and is due to finite-size effects in small systems. 
From the intersections of sizes L > 8, a value for the crit- 
ical randomness of /ij?'"'^"'' = 4.18(2) is estimated. Due to 
the finite-size effect observed here, small system sizes will 
be excluded from subsequents fits if they don't match the 
leading behavior. 

To observe the specific heat singularity, here the bond 
energy Ej is studiedt^, given by 



E, 



(5) 



where the sum is over nearest- neighbor pairs. Ej has an 
energy-like singularity in the vicinity of the phase bound- 
ary. Having differentiated analytically with respect to J, 
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now J = 1 is set, and T = considered only. A spe- 
cific heat-like quantity is obtained by differentiating Ej 
numerically with respect to the random field h. A first- 
order finite difference is used to determine the derivative 
of Ej numerically and, since this is a more accurate repre- 
sentation of the derivative at the midpoint of the interval 
than at either endpoint, the "specific heat", C, at T = 
is defined to be 

ffn + h2 \ _ [Ej{h,)U-[Ej{h2)U ... 

where hi and /12 are two "close-by" values of h. A suf- 
ficiently fine mesh of random-field values is chosen such 
that the resulting data for C is smooth. Error bars are 
obtained by determining the "specific heat" from the cor- 
responding finite difference as in Eq. (||) for each sample 
separately, and computing the standard deviation. The 
error bar is, as usual, the standard deviation divided by 

\/A^samp 1- 



3 



2 



1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 


1 




« 4 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




G— 

■ 

_ A- 
X- 

- *- 


-0 L=32 §f 
■ a L=24 f i 
^L=16 IT ^ 
-a.L=12 
-«L=8 b I 
-*L=4 «( W2 


- ^^'\ - 

-3'^% ■ 






4.1 4.2 43 44 




J/ \ 










\ 

**^ 




_1 1 L. 


J 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 i_ 


J l_ 


3 


4 


5 


f 



h 

FIG. 4: Specific lieat-like quantity C = d^FjdhdJ as a func- 
tion of the randomness h for L = 4, 8, 16, 32. Error bars are 
only shown for L = 32. For the other system sizes the error 
bars are even smaller. The inset shows the region near the 
peaks enlarged for L > 8. Lines are guides to the eyes only. 

In Fig. 1^ the numerical results are exposed. The "spe- 
cific heat" exhibits clear peaks, which grow in system size 
and move slightly left. The number of samples used is 
larger near the peak to compensate for the greater sample 
to sample fluctuations in this region. 

In a finite system, finite-size scaling predicts for the 
singular part Cg of the specific heat 

~L"/'^c((/i-/ie)Li/'') , (7) 

where v is the correlation length exponent. The "specific 
heat" peak will occur when the argument of the scaling 
function C takes some value, a\ say, so the peak position 
h*{L) varies as 

h*{L)-h,^aiL-^'\ (8) 



and the value of the singular part of the "specific heat" 
at the peak varies as 

C^"i'^'^(L) - L"'". (9) 
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FIG. 5: A plot of the random field where the "specific heat" 
attains its maximum, as a function of system size L. The 
solid line shows a fit to the function h*{L) = he + a\L~^^" 
with he = 4.18, 1/u = 0.78, and ai = 1.34. The inset shows 
the data as a function of L^^^" . 

For each system size, parabolic fits were performed in 
the region of the peak to obtain h* (L) and the height of 
the peak, ^^^^(L). The shift of the maximum according 
to Eq. (j^) can be used to estimate the infinite-size criti- 
cal strength of the random field, he and the correlation- 
length exponent v. The best fit gives 

= 4.182 ±0.006, 1/1/= 1.28 ±0.14, (10) 

see Fig. |. The probability Q was determined that 
the value of = X^^i C^'"/.^'^''' with N data points 
(xi, Ui ± (Ti)-|fitted to the function /, is worse than in the 
current fitc3 to quantify the quality of the fit. Here a 
value of Q = 0.06 was obtained , which is not very good. 
The reason is, that the error bars used for the positions 
and heights of the maxima were only the statistical error 
bars obtained from the fit of the parabolas, and which 
are often surprisingly small. Systematic errors, resulting 
from the fact that the peaks are in fact not parabolic, 
are not included in this way. Therefore, fits to a fourth 
order polynomials were tried, but the results turned out 
to be very unreliable and the fits very unstable against 
the change of the window over which the data was fitted. 
Hence, the parabolic fits were kept, were these effects 
were smaller, and the final value quoted is he = 4.18(1). 

Next, the singular behavior of C is determined by an- 
alyziag, how the peak value C'™^'^ scales with L. Please 
noteta that at T = the singular behavior C is equal 
to the singular behavior of C = ~d^F/dh^, since C = 
-C'h/J at r = 0. 
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FIG. 6; The maximum (7™^" of the "specific heat" as a func- 
tion of system size with logarithmically scaled L-axis. The 
lines shows the function a + blogL with a = —0.53 and b = 
1.33. The inset shows the data in a double-logarithmic plot, 
the solid line being the function aiL°'^'^ with a/u = 0.419 and 
a2 = 0.98. 

The first hypothesis is that the "specific heat" di- 
verges logarithmically, as found in experimentsJor three- 
dimensional diluted antiferromagnets in a fieldE3. Hence, 
the values of C^'^^ are fitted {L > 8) to a function of the 
form 

C'^'^^[L)=a + b\ogL, (11) 

where the constant term a comes partly from the regular 
piece of the "specific heat" . The quality Q = 0.23 of the 
fit is fair. In Fig. |^ the "specific heat" is shown in a 
logarithmic plot together with the fit function (pH). For 
larger system sizes the datapoints follow very nicely a 
straight line, suggesting that C indeed may diverge log- 
arithmically. Please note that the "specific heat" peaks 
seem to be rather symmetrical, which means that the 
amplitudes A^, (for h > hch < he) oi C are almost 
equal. This is exactly jjihat one expects in the case of a 
logarithmic divergenceEil. 

Also the possibility of a algebraic divergence was 
tested. For this purpose the data was fitted (i > 8) 
to the function 

C'°^'^(L) = aai"/'', (12) 

resulting in a/v = 0.419(9) and 02 = 0.98(2), the result 
is shown together with the data in a double logarithmic 
plot in the inset of Fig. ||. The quality Q = 9 x lO"'^ 
of the fit is very bad. To check whether this is an ef- 
fect of including too small system sizes into the fit, also 
a fit using only system sizes L > 12 was performed, re- 
sulting in a/v = 0.416(6), 02 = 0.99(1) and Q = 0.16 
which is much better. Since for i > 12 the logarith- 
mic fit has also a better quality Q = 0.55, and because 
the negative curvature in the double-logarithmic plot is 



more pronounced than a possible positive curvature in 
the single-logarithmic plot, still a logarithmic divergence 
seems more likely from this data, i.e. 

a = 0, (13) 

but an algebraic behavior with a small exponent cannot 
be excluded. 

Next, the critical behavior of the magnetization is stud- 
ied. The predictions from finite-size scaling is that near 
the critical point 

m{h) = L-^'^miih - he)L^'''). (14) 

This means, that by plotting m{h)L^/'^ against L^/'^{h — 
he) with correct parameters hc,v and P/v, the data 
points for different system sizes should collapse onto a 
single curve near {h — he) ~ 0. The values he — 4.18 and 
1/v = 1.28 from above were used. With 

= 0.17(5) (15) 

the best collapse of the data for L > 8 was obtained, 
which is presented in Fig. |^. 




FIG. 7: Scaling plot of the rescaled absolute value m{h)L^^" 
magnetization as a function of of {h — hc)L}^'' with he = 4.18, 
1/v = 1.28 and P/u — 0.17. Error bars are smaller than 
symbol sizes. Lines are guides to the eyes only. 

The singular behavior of the disconnected susceptibil- 
ity 

can be obtained in an analogous way to the magnetiza- 
tion. The following scaling behavior is assumed: 

XdUh) = i"/"Xdis((/i - he)L^/n- (17) 

From collapsing the date curves for L > 8, using he = 
4.18 and l/i^ = 1.28, a value of 

7/1/ = 3.63(0.05) (18) 
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was found. The scaling plot is shown in Fig. g. Please 
note that the scaling behavior oih < is worse than for 
h > he (also, to a lesser extent, in Fig. 0). The reason is 
that for smaller fields the systems quickly become fully 
ordered {m = 1), i.e. scaling does not hold. 
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FIG. 9: Susceptibility x a^s a function of the random-field 
strength h for system sizes L = 4, 6, 8, 12, 16, and 24. Error 
bars (shown for L = 24) are smaller than symbol sizes. Lines 
are guide to the eyes only. 



FIG. 8: Scaling plot of the rescaled disconnected susceptibility 
XdiB^L'^" as a function of of {h - with = 4.18, 

l/v — 1.28 and '^/v = 3.63. Error bars are smaller than 
symbol sizes. Lines are guides to the eyes only. 

Finally, the susceptibility and its related critical expo- 
nent 7 is determined. This is done by considering the 
response to a small uniform external field H > 0. For 
each realization, the sign of H is chosen in the direction 
of the magnetization of the ground state. This prevents 
the whole system from flipping when applying a mag- 
netic field to a system which is almost ferromagnetically 
ordered, which would mimic a false high susceptibility. 
The scaling behavior of the magnetization should not be 
affected by this choice. 

The ground state for each realization and each value of 
h is calculated for H„ = 0,Hl,2HlAHl- Near H = 0, 
the data points can be fitted very well with a parabola, 
the coefficient of the linear term gives the zero field sus- 
ceptibility 

X - dm/dH\H=a. (19) 

To cope for the expected strong increase of x with the 
system size, must strongly decrease with L, in the 
order of the expected increase. For details see Ref. |l^, 
the values of are shown in Tab. ||. Then, for each sys- 
tem size, a fit to a parabola through the four data points 
for the average magnetization to(77„) isi^erformed. To 
estimate the error, a jackknife analysis^ was used, in 
which the results for the magnetizations (for each sys- 
tem size and each strength of the disorder) is divided 
into K blocks, the average values calculated K times, 
each time omitting one of the blocks, and then K fits are 
performed. The error bar is estimated from the variance 
of the K results for the linear fitting parameter. Here 



K = 50 was used and checked that the result does not 
depend much on the choice of K. 

The susceptibility x as a function of h is presented 
in Fig. ^ for selected system sizes. It is seen that the 
height of the peak grows much faster than for the "spe- 
cific heat". For the susceptibility, the following scaling 
behavior is expected: 

x{h)^L'''^^{{h-h,)L^/^). (20) 

To analyze the divergence of x, again parabolas were fit- 
ted to the data points near the peak to obtain the posi- 
tions h*{L) and x"^^^{L) of the maximum. By fitting the 
data for i > 8 to a function x'°''''(i) = asL^/", where 
7 describes the decay of the "connected" correlations at 
criticality, the following values were obtained {Q — 0.54) 

7/z^ = 1.82 ±0.01, (21) 

see Fig. ^ 

It was also tried to fit the positions H*{L) of the max- 
ima of the susceptibility to a scaling function of the from 
(^). But the quality of the fit was very bad {Q < 10"'') 
for all ranges [Lmin, 24] considered and the result for he 
was always too large {he > 4.26 ±0.03). This value is not 
only not compatible with the result from the "specific 
heat" C, but in clear contradiction to the Binder param- 
eter (see Fig. ||). The reason is probably the stronger 
finite-size dependence of the position of the susceptibil- 
ity peak in comparison to the position of the "specific 
heat" peak, please compare Figs. ^ and ^. Also, due 
to the four times higher numerical effort to determine x, 
only simulations for L < 24 were performed. Further- 
more, the peaks for the susceptibility are much broader 
than for the "specific heat" , hence it is harder to deter- 
mine the position of the peak precisely. Thus, the result 
he = 4.18(1) from the previous measurement was kept. 
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FIG. 10: The maximum i^™^" of the susceptibility as a func- 
tion of system size L in a double logarithmic plot. The solid 
line represents a fit to the function ^'^^{L) — a^l?^^ , for 
sizes i > 8 yielding -^jv = 1.82 and az — 0.066. 

IV. DISCUSSION 

Using graph-theoretical algorithms, exact ground 
states of random-field Ising systems were calculated in 
polynomial time. Using the LEDA library, system sizes 
up to = 32'* could be considered. 

All critical exponents, describing the order-to-disorder 
transition at he were calculated independently. The fol- 
lowing values, using 7 = i/(2 — 77) and 7 = !^(4 — rj), for 
the values of the critical field he and the exponents were 
obtained: 

/ic = 4.18±0.01, = 0.78 ±0.10 

a ^ /3 = 0.13 ±0.05 ,^2) 

T] =0.18±0.01 ry= 0.37 ±0.05 

7 =1.42±0.20 7 = 2.83 ±0.50. 

Please note that the values for (3, 7 and 7 carry a large 
error bar due to the uncertainty in the correlation-length 
exponent v. 

The results for the critical field, the correlation expo- 
nent and the exponent of the magnetization are com- 
patible with values found formerly-via the exact ground 
state calculations of small systema23 up to L = 10, where 
he = 4.17(5), P = 0.13(2) and ly = 0.8(1) were obtained. 
The result for the susceptibility exponent is compatible 
with the result 7 = 1.45.(5) which was found in a high- 
temperature expansionll3. The results obtained here axa 
also compatible with the exponents obtained recentlj^ 
by the use of exact ground states as well but by evaluating 
mainly other quantities like distributions of domain-wall 
energies or fractal properties of domain walls. The main 
difference is that the results from Ref. |3^ have a slight 



preference for the exponent a of the "specific heat" to be 
positive but small. 

Next, the validity of the scaling relation for the "spe- 
cific heat" is checked. For the Rushbrooke equality 

a + 2/3 + 7 = 2 (23) 

one gets a + 2/3 + 7 = 1.68(30), which is not very good 
but still within the error bars almost fulfilled. In case 
the algebraic divergence is taken (a = 0.33(5)) a value 
of a + 2/3 + 7 = 2.11(35) is obtained, which fulfills the 
equation better. 

The deviation of the hyper-scaling relations from the 
pure case is obtained by replacing d hy d — 9, with the 
exponent = 2 — rj+rj = 1.81(6). E.g. the hyper-scaling 
relation 

2~a^v{d-e) (24) 

is also fulfilled within error bars {v{d — 6) = 1.70(30)), 
again the case were it is assumed that the "specific heat" 
diverges faster than logarithmic matches the relation bet- 
ter. 

Finally, we turn to the question whether there are two 
or three independent exponents. In the case of two ex- 
ponents, the Schwartz-Soffer equationliJ 

7 = 27 (25) 

holds, which is compatible with the result found here. 
Hence, the two-exponent scenario is supported. 

To summarize, all critical exponents of the four- 
dimensional RFIM were determined independently. All 
"classical" (hyper-) scaling relations are fulfilled and the 
two-independent-exponents scenario is supported. The 
scaling relation proposed in Ref. |2^ seems not to be ful- 
filled. The largest uncertainty in the results presented 
here is in the value of a, whether a = or a > and 
small. To finally decide this question, much larger system 
sizes must be studied, which are currently out of reach. 
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